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We study relaxation and dephasing in a strongly driven two-level system interacting with its 
environment. We develop a theory which gives a straightforward physical picture of the complex 
dynamics of the system in terms of dressed states. In addition to the dressing of the energy diagram, 
we describe the dressing of relaxation and dephasing. We find a good quantitative agreement between 
the theoretical calculations and measurements of a superconducting qubit driven by an intense 
microwave field. The competition of various processes leads to a rich structure in the observed 
behavior, including signatures of population inversion. 
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An essential question in quantum theory is how a sys- 
tem is affected by its interaction with its environment. 
There has been great progress in recent years describing 
this interaction through decoherence theory, which quan- 
tifies the effects in terms of relaxation and dephasing. 
An important question in this field is how the interac- 
tion of the system with its environment is modified when 
the system is strongly driven. There has been a signifi- 
cant theoretical effort to understand this problem pQ , but 
there remains a variety of different theoretical approaches 
with few experimental results to guide progress. This 
has changed recently with the emergence of the new field 
of circuit quantum electrodynamics [2j |3] , where nano- 
electronic circuits interact with photons at the quantum 
level. The design flexibility afforded by these solid-state 
systems has allowed the exploration of new regimes of 
drive strength and new forms of interaction [4j [5j |6] . In 
this Letter, we use the techniques of circuit QED to make 
a quantitative comparison between the theory and exper- 
iment of relaxation and dephasing in a strongly driven 
system. 

We present a microscopic model of how the driven 
system interacts with the quantum vacuum of the en- 
vironment. We exploit a hierarchy of energy scales to 
develop an analytic description which gives a straight- 
forward physical picture of this complex system in terms 
of dressed states. We show that by including a mini- 
mal number of parameters, which describe the spectral 
density of the vacuum, we can explain a wide variety 
of observed effects, including population inversion in the 
dressed states. 

Dressed states of superconducting circuits have re- 
cently received attention in the context of quantum infor- 
mation 0. In particular, theoretical work on quantum 
state detection, i.e., qubit readout [8 , has suggested that 
relaxation and dephasing of the dressed states may be an 
important factor limiting performance [9]. 

Our artificial atom, the single Cooper-pair box (SCB), 
is composed of a superconducting Al island connected 
to a superconducting reservoir by a small area Joseph- 
son junction [10 j. The two charge-basis states of the 



SCB represents the presence (absence) of a single ex- 
tra Cooper-pair on the island. The Hamiltonian of the 
SCB coupled to the driving microwave field is H = 
-\E C h°z - \Ejcf x + hw^a + ga z (a + a 1 "), where cr* 
are the Pauli spin matrices and a\ a are the creation 
and annihilation operators for the microwave field. The 
first two terms represent the uncoupled SCB Hamilto- 
nian, where Ecu — Eq(1 — 2n g ) is the electrostatic en- 
ergy difference between the ground and excited state of 
the qubit and Ej is the Josephson coupling energy. Here 
E Q = (2e) 2 /2C E is the Cooper -pair charging energy, 
is the total capacitance of the island, and n g = C g V g /2e 
is the dc gate charge used to tune the SCB. We contact 
the island with a two junction SQUID which allows us 
to tune Ej with a small magnetic field. The third term 
represents the free driving field, and the last term rep- 
resents the coupling, with strength g, between the SCB 
and the microwave amplitude. 

We measure the dressed states by coupling the driven 
SCB to an rf oscillator (Fig. []Je)). We probe the oscil- 
lator with a weak rf field, measuring the magnitude and 
phase of the the rf reflection coefflent, Sn. In a typical 
measurement, we use an external magnetic field to fix 
a value for Ej. We then produce a 2D map of Sn by 
slowly sweeping n g while stepping the microwave ampli- 
tude. For more details of the experimental setup, see [4]. 
In Fig. [TJa) & (b), we present measurements of Sn for 
representative values of Ej/h — 2.6 GHz, Eq/H = 62 
GHz, w m /2tt = 7 GHz and c^ r/ /27r = 0.65 GHz. We see 
a rich response in both the magnitude and phase. 

To understand the data, we start by considering the 
Hamiltonian of our system, H, in a region of n g where 
Ech ~ nhuj^. To first order, we get a ladder of energies 

£$ 0) = Nhu>» ± \y](E C h-nhwtf + (A„(a)) 2 

with pairs of an "excited" state, |l(iV)), and a "ground" 
state, |0(iV)), that repeat for different photon numbers 
N (Fig. [iff)) QH H2II3|. Here A n (a) = EjJ n (a) is 
the dressed gap between these states, which varies with 
the normalized microwave amplitude a = 2EQn /Jj /huj fI , 
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FIG. 1: (Color Online) RF response of the strongly driven SCB coupled to an oscillator. The images represent the rf reflection 
coefficient of the system, 6*11, as a function of dc gate charge, n g , and microwave amplitude, A^. (a) and (b) show the data 
for the phase, arg(5n), and magnitude, \Sn\ — 1, respectively, (c) and (d) show the corresponding theoretical calculations. 
The fitting parameters represent a model for the environmental charge noise, sketched in (g), which is responsible for dressed 
relaxation and dephasing. The data and theory are plotted with the same color scales, (e) Schematic of the SCB and resonator, 
(f) Energy level diagram showing the dressed levels and transitions, (g) Sketch of the environmental spectral density showing 
the 1/ f tail, oscillator peak, and Ohmic background. 



n M = ^^A^jle being the microwave amplitude in units of 
2e. Here is the microwave amplitude at the generator 
and 7^ is the microwave coupling. Finally, J n is a Bessel 
function of order n. 

We have previously shown that the absorption features 
in Fig. [IJb) arise from the resonant interaction of the 
dressed states and the readout oscillator [4 j. However, 
the phase response is more varied. Referring to Fig. 
[TJa), we see that there is a cone at high power where 
the response is relatively simple, showing vertical stripes 
where the phase shift is unimodal. We showed this re- 
sponse could be explained well by the dispersive shift of 
the oscillator frequency caused by its coupling to the near 
resonant dressed states, an effect related to the quantum 
capacitance of the states p~4j [15]. At lower powers, we 
see that there are regions where the phase shift becomes 
bimodal, changing signs as a function n g when moving 
across a resonance. This is seen clearly in Fig[2|a) as a 
change in color from red to blue. We note that the tran- 



sition in the character of the response does not occur at 
a simple uniform threshold. 

To explain this rich variety of phenomena, we must un- 
derstand how relaxation and dephasing are dressed in our 
strongly driven system. We start by considering transi- 
tions between states in this ladder which are induced by 
the system's coupling to the environmental bath, seen 
through the oscillator. We do this by writing the master 
equation describing the dynamics of the density matrix, 
p, which spans our set of states. The master equation 
contains a set of transitions rates which are proportional 
to the square modulus of the matrix elements of the oper- 
ator which couples the system to the bath. To calculate 
these matrix elements, we first extend the calculation of 
our dressed states to next order in Ej/mhw^. We in- 
clude the contribution of all off-resonant states at this 
order, including the effects of m-photon transitions for 
arbitrary values of m. After calculating the rates, we 
then reduce the full master equation to a master equa- 
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tion for an effective two- level system by tracing over N. 
It is this reduced master equation that we use to describe 
the dynamics of the reduced dressed states, |1) and |0), 
interacting with our readout oscillator. 

We first consider the effects of charge noise in the 
environment, which is typically the dominant source of 
noise for a SCB. The noise couples through a z and has 
a spectral density Sq(lu). We start by assuming that 
the bath is at zero temperature. The reduced mas- 
ter equation then includes three rates: the relaxation 
Z)m>o-^rez5 tne excitation rate, 
; and the dephasing rate, 1^(77) = 
(p 1 ^m>() L ^ • Here r] is the mixing angle across the 
dressed degeneracy point defined by tan r] = A n /(Ech — 
mhujmu). The different rates are TfJ = A 2 sin 2 tjSq(A) 
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where A = 1 - (E J /2hu; fI ) 2 ^ m _, J m-n(«)M 

= E 2 SQ(mhx)^) I ' (mhw jJi ) 2 . We can include the ef- 
fects of a finite bath temperature, T, by adding to r ea;c 
the term T re i exp(— ftA/fe^T). A few comments are in 
order. First, we note that T^J and Y s * d are the standard 
rates we would calculate for an undriven two-level system 
coupled to a bath (assuming A ~ 1). These rates repre- 
sent transitions that do not change N. The other rates, 
r m , represent transition that change N by m photons. 
Next, we note that we have an effective excitation rate 
even though we have only included the effects of charge 
relaxation in the full master equation. This is related 
to transitions that take us from a ground dressed state, 
|0(iV)), to an excited dressed state with fewer photons, 
\1(N — m)) (FigjTJf)). In this transition, the total sys- 
tem loses energy to the bath, but our effective two-level 
system appears to be excited. 

We model the readout of the dressed states in the 
following way. Starting from the reduced master equa- 
tion, we derive the semiclassical Bloch equations for the 
dressed state charge driven by the small rf probe volt- 
age V r f = 2en r f /C r f. Besides the three rates above, 
the Bloch equations also include an effective tempera- 
ture, T e ff, which determines the equilibrium occupation 
of the states. This is calculated as a function of a and rj 
as 



T eff = 
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We then take the standard solution of the Bloch equa- 
tions for the in-phase and quadrature component of the 
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FIG. 2: (Color Online) (a) Dressed bloch response of the 1- 
photon dressed state. The (color) scale is the same as Fig. [l] 
The two left panels are the theory and data for higher am- 
plitude, while the right panels are for lower amplitude, (b) 
The extracted values of F^J as a function of the detuning 
between the dressed states and oscillator, for positive detun- 
ing, (c) The calculated values of T exc (r] = tt) and T re i(ri — tt) 
as a function of the microwave amplitude A^ . 



(excluding FfJ) 



dressed state charge, Qi and Q q . These are used to cal- 
culate an effective parallel resistance and capacitance, 
Reff = V rf /(3Q q = V rf /2nf c /3Q g and C eff = /3Qi/V rf , 
which are then used to calculate the reflection coefficient 
Sii. This is essentially the same procedure used in [4] 
except that we now include the calculated T e f f . In addi- 
tion, the coupling constant /3 = C r f / has been added 
to better account for the coupling of the dressed state 
charge to the oscillator. 

To compare our theory with experiment, we need to 
have a model for the environmental spectral density 
Sq{uj) (see Fig jib). The typical starting point is to 
assume that Sq(u) is Ohmic with an additional con- 
tribution from 1/f charge noise [17]. Also, because of 
the readout oscillator, Sq(co) will not be smooth around 
uj r f. For the purposes of fitting, we will then describe 
the environment by three parameters: Sq(uj w 0) which 
determines r^ £d , Sq(uo = A n ) which determines r**f, 
and a high-frequency coupling constant ft, such that 
Sq(lj » A n ) oc k 2 cuRo. k then determines the rates for 
the m-photon relaxation processes. We note that com- 
pared to the simple model in [4] , we have only added one 
fit parameter, which is Hi. 

In Fig. |2|a), we show the results of performing a de- 
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tailed fit to the one-photon resonance of the data in Fig. 
[I] The magnitude and phase data are fit simultaneously, 
which is important for the parameters to be constrained, 
although we only show the phase data. For compact- 
ness, the vertical axis is split and the higher and lower 
lobes are plotted side- by-side. We see that the agreement 
between data and theory is very good for both lobes of 
the response, despite the fact that the experimental re- 
sponses look very different. The three parameters men- 
tioned above are allowed to vary independently for each 
value of a. However, we find that the variation of the 
value of k is less than 10%, which is comparable to the 
random error in each fit. The extracted value of Sq(Huj I j j ) 
would translate into a relaxation time for a typical charge 
qubit of 1/TreZ ~ 300 ns, consistent with observed values. 
In Fig. p[b), we plot the extracted values of the r£*f, 
proportional to 5g(A n ). We see that Sq(co) is peaked 
around the oscillator frequency uj r f, as we expect and as 
we observed before [4 . It is worth noting again, that de- 
spite the very different appearance of the response in the 
two lobes, the extracted rates are consistent with each 
other, falling on the same curve. Taken together, these 
results confirm that we have a good understanding of how 
charge relaxation takes place in the extended ladder of 
dressed states, and how that reduces to relaxation and 
excitation in the reduced basis |1) and |0). 

It is possible to understand the response in straightfor- 
ward physical terms within our dressed state interpreta- 
tion. First of all, we can understand the bimodal nature 
of the phase response at lower a (right side of Fig. [2|a)). 
If we consider the n-photon resonance on one side of the 
dressed degeneracy, near 77 ~ 0, the dominate terms in 
the rates (eqs. [l][2| are T re i(rf) oc Jq (a) ~ (1 — a) 2 and 
r eX c(v) oc Jl +1 {a) ~ a 2(n+1) . Clearly then, for a « 1 
relaxation dominates. On the other side, at 77 ~ 7r, we 
find instead T exc {rf) oc J^(a) and Y re i(r]) oc J^ +1 (a), im- 
plying that excitation is in fact dominant. This implies 
that the reduced states |1) and |0) will become inverted in 
this region. At the degeneracy point, T re i(rj) ~ F exc (r]) 
(if we ignore the standard relaxation for now) and we 
expect the populations to equalize. The Bloch equations 
tell us that the phase response of the excited state has 
the opposite sign from the ground state response. There- 
fore, as we move across the dressed degeneracy point, we 
expect the phase shift to start with one value, go to zero, 
and then change sign. This is exactly the bimodal shape 
that we observe. 

As we move to higher values of a, the situation 
changes. The dominant Jo term in the rates decays 
sharply, leaving the summed contributions from many 
different photon transitions. As shown in Fig. [2|c), these 
summed rates become essentially equal for larger ampli- 
tudes. Essentially, the strongly driven transition satu- 
rates. Once this happens, the dynamics of the reduced 
dressed states are determined only by T^J. This is why 
we see a crossover to the simple form of the response at 



higher a. 

Having shown that we can fit the response for partic- 
ular values of a, we can also wonder if it is possible to 
explain the response over the whole range of parameters. 
To do this, we need to account for additional dephasing 
which is not captured by (|3|. This is not in fact surpris- 
ing since the first order contribution should vanish at the 
dressed degeneracy points. We use an adiabatic approx- 
imation to model the second order charge noise through 
a pure dephasing term connected to an effective <j x oper- 
ator [18 . To lowest order in Ej/mftw^, we find that the 
contribution of m-photon transitions is 

r™ x (v) = sin2 V {j-{m+n) + Jm-n) 2 S X {mhw^) (4) 

where Sx(w) is then the spectral density of the noise 
associated with the effective a x operator. 

In Fig. [TJc) & (d), we show calculations of the re- 
sponse covering the full range of data. We no longer fit 
the parameters associated with relaxation in these calcu- 
lations. Instead, we model the environment based on the 
parameters extracted from the fits presented in Fig. [2] 
In modeling the additional dephasing, we have assumed 
that the spectrum has an effective cutoff above fio^, i.e. 
we only include the m = 0, 1 transitions. This is justified 
by the adiabatic approximation mentioned above. There- 
fore, there is a total of just 2 fitting parameters for both 
2D plots combined, Sx(0) and Sx{fouy)- We see that the 
agreement for the magnitude data is striking and that it 
is also quite good for the phase data. In particular, we 
reproduce both the quantitative level of the response and 
all the qualitative features mentioned earlier. (The pos- 
sible exception to this statement is the 0-photon phase 
response at small A^. This is not surprising, however, 
as in this region the dressed gap is far detuned from the 
resonator and our model for the response does not ap- 
ply.) Taken as a whole, this is strong evidence for our 
interpretation. 
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